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Abstract 

Large-scale peculiar motions are believed to reflect the local inhomogeneity and anisotropy 
of the universe, triggered by the ongoing process of structure formation. As a result, realistic 
observers do not follow the smooth Hubble flow but have a peculiar, ‘tilt’, velocity relative 
to it. Our Local Group of galaxies, in particular, moves with respect to the universal expan¬ 
sion at a speed of roughly 600 km/sec. Relative motion effects are known to interfere with 
the observations and their interpretation. The strong dipolar anisotropy seen in the Gosmic 
Microwave Background, for example, is not treated as a sign of real universal anisotropy, but 
as a mere artifact of our peculiar motion relative to the Hubble flow. With these in mind, 
we look into the implications of large-scale bulk motions for the kinematics of their associ¬ 
ated observers, by adopting a ‘tilted’ Friedmann model. Our aim is to examine whether the 
deceleration parameter measured in the rest-frame of the bulk flow can differ from that of 
the actual universe due to relative-motion effects alone. We find that there is a difference, 
which depends on the speed as well as the scale of the bulk motion. The faster and the 
smaller the drifting domain, the larger the difference. In principle, this allows relatively slow 
peculiar velocities to have a disproportionately strong effect on the value of the deceleration 
parameter measured by observers within bulk flows of, say, few hundred megaparsecs. In 
fact, under certain circumstances, it is even possible to change the sign of the deceleration 
parameter. It goes without saying that all these effects vanish identically in the Hubble 
frame, which makes then an illusion and mere artifact of the observers’ relative motion. 


1 Introduction 

On large enough scales, of the order of several hundred megaparsecs and beyond, our universe 
appears uniform to a very good approximation. On sufficiently small scales, however, there is a 
lot of structure (in the form of galaxies, galaxy clusters, superclusters and voids), which destroys 
local homogeneity and isotropy to a larger or lesser degree. Bulk peculiar motions are a direct 
result of the ongoing process of structure formation. Typical galaxies do not follow the smooth 
Hubble expansion but have a ‘drift’ motion relative to it [T]. Our Local Group of galaxies, for 
example, has a peculiar velocity of approximately 600 km/sec with respect to the frame of the 
Cosmic Microwave Background (CMB). The latter is defined as the coordinate system where 
the CMB dipole vanishes and provides the reference frame of the universal expansion, relative 
to which we can define and measure large-scale bulk motions. 
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Peculiar velocities are expected to fade away as we move on to progressively larger wave¬ 
lengths. The current concordance cosmological model, namely the ACDM paradigm, allows for 
peculiar velocities up to 250 km/sec on scales close to 100/h Mpc, with decreasing magnitude as 
we move on to larger lengths 0 Taking cosmic variance into account could bring the aforemen¬ 
tioned limit higher, though not to the magnitudes reported by a number of recent surveys. The 
work of [2], in particular, claims bulk velocities about 400 km/sec within a sphere of lOO/Zi Mpc 
size. Using different methods, however, other authors have reported lower amplitudes for the 
local bulk flow laai, some of which lie within the general ACDM predictions [5]. Much faster 
peculiar motions, between 1000 km/sec and 4000 km/sec, covering scales of Gpc-order have also 
been reported in the literature [6]. Nevertheless, such very large-scale ‘dark flows’, as they are 
commonly known, appear to be at odds with the Planck constraints. The latter set an upper 
limit of roughly 400 km/sec within regions larger than 500 Mpc [9], but the issue may not be 
over yet m- 

Relative motions are known to interfere with the observations and with their interpretation. 
The strong dipolar anisotropy seen in the CMB spectrum, for instance, is not treated as a sign of 
real cosmic anisotropy, but as a mere artifact (a Doppler-like effect) of our motion with respect 
to the smooth Hubble expansion. With these in mind, we will employ a tilted Friedmann- 
Robertson-Walker (FRW) cosmology and a two-fluid approach to investigate the implications of 
large-scale bulk motions for the mean kinematics of the ‘drifting’ observers. We introduce two 
timelike 4-velocity fields, tangent to the worldlines of their corresponding observers, which move 
relative to each other. The first 4-velocity held will be identihed with the smooth Hubble how 
and will dehne our frame of reference. Relative to this coordinate system, we will assume that the 
universe is an FRW spacetime hlled with conventional pressureless matter. The second 4-velocity 
held, which has a small but hnite peculiar velocity with respect to the hrst, will dehne the rest- 
frame of a realistic observer in a typical galaxy that moves relative to the Hubble expansion. 
Treating the peculiar velocity as a small fraction of the Hubble speed, on the corresponding 
scale, conhnes the domain of our study to wavelengths larger than ~ 100 Mpc. The latter is also 
thought to mark the ‘homogeneity threshold’ of our cosmos, namely the scale beyond which the 
universe behaves as an almost-FRW model, although this is a working hypothesis rather than a 
hard observational fact. 

Assuming that the drifting observers reside within a large-scale bulk flow, like those reported 
in m-m, we focus on their average volume expansion and more specifically on the deceleration 
parameter measured in their own rest-frame. Following [ZllH], we find that bulk motions can 
modify the local expansion (Hubble) rate and that the magnitude of the change depends on 
the speed of the peculiar flow. Not surprisingly, and in agreement with 13 i, we also find 
that variations in the local volume expansion change the associated deceleration/acceleration 
rate and therefore the local value of the deceleration parameter. To estimate these changes, we 
extend the work of 1318] and turn to linear cosmological perturbation theory, which allows us to 
connect peculiar velocities to density perturbations. The result is a simple relation expressing 
the deceleration parameter, as measured in the rest-frame of an observer living inside a large- 
scale bulk flow, in terms of the deceleration parameter of the actual universe and gradients of 
the peculiar velocity field. 

^The dimensionless variable h — H/IQQ km/sec Mpc defines the normalised Hubble parameter 


2 



Perhaps the most interesting feature of the aforementioned mathematical expression is its 
scale-dependence. We find that the effect of relative motion on the deceleration parameter, as 
measured in the rest-frame of the bulk flow, increases with decreasing scale. The threshold is the 
Hubble radius, which makes coherent drift motions with sizes of few hundred Mpc particularly 
‘vulnerable’ to relative-motion effects. On these scales, even the peculiar velocities allowed by 
the ACDM model can seriously affect the local value of the deceleration parameter. In addition 
to the magnitude of the velocity and the size of the bulk motion, the impact on the deceleration 
parameter also depends on whether the peculiar flow is (slightly) contracting or expanding. 
Crucially, when the volume of the bulk is decreasing, relative motion effects can even change 
the sign of the deceleration parameter there. None of these effects is ‘real’, of course, since they 
all vanish in the CMB-frame. They are mere artifacts of the observers’ peculiar motion, just 
like the CMB dipole, but ignoring them could potentially lead to a serious misinterpretation of 
the observational data. Put another way, measuring a negative deceleration parameter in the 
rest-frame of a bulk flow may not necessarily mean that the universe itself is accelerating. 

2 Aspects of relative motion 

Consider two families of relatively moving observers, with worldlines tangent to the (timelike) 
4-velocity fields Ua and Ua respectively. Assuming that Va is the peculiar velocity of the latter 
family with respect to the former, we have 

Ua = l{Ua + Va) , (1) 

where 7 = (1 — > 1 is the associated Lorentz factor. Note that UaV°‘ = 0, 

and UaU^- = — 1 = UqU® by construction (e.g. see m)- The above relation also determines the 
‘tilt’ angle between the two 4-velocity fields. In particular, the inner product UaU°‘ = —7 < —1, 
defines the hyperbolic angle (/3) between Ua and Ua (see Fig. d]), so that cosh /3 = —UaU°' = 7 > 1 
and /3 = ln( 7 -|- 1 / 7 ^ — I). In what follows we will only look at non-relativistic peculiar motions, 
with <C 1 and 7 ~ 1. Therefore, in our case, the Lorentz boost reduces to 

Ua — Ua -|- Va , ( 2 ) 

while the tilt angle between Ua and Ua becomes very small (i.e. /3 ~ In 1 ~ 0 ). 

The velocity fields Ua and Ua define the time-direction of their associated observers, while 
their 3-dimensional rest-spaces are normal to the corresponding 4-velocities. Spatial geometry 
is determined via the so-called projection tensors. These are defined by hab = gab + UaUb and 
hab = 9ab + UaUb, where Qab is the metric of the whole spacetime, so that habU^ = 0 = In 

addition, hab = ^(afe); ^ab = hach^b and /ia“ = 3, with exactly analogous relations holding in the 
‘tilted’ frame as well (e.g. see [Kills]). Finally, employing the aforementioned 4-velocity fields 
and their corresponding projectors, we can define temporal and spatial derivatives, relative to the 
two frames. Thus, from now on, overdots and primes will indicate time differentiation in the Ua- 
frame and its tilted counterpart respectively (i.e. ' = u“Va and ' = u°'S/a, with Va representing 

^Irrotational 4-velocity fields are hypersurface orthogonal as well. Then, the projector also acts as the metric 
tensor of the observers’ 3-D space. In this study the Wa-field will be irrotational by construction (see § [4] below). 
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Figure 1: Observers with 4-velocity Ua have peculiar velocity Va = vca relative to the Ua- 
field. By construction, Uae“ = 0 = and eae“ = 1 = eae°‘. The hyperbolic angle /3, with 
cosh/3 = —UaU^, defines the ‘tilt’ between Ua and Ua, while S is the 3-D space normal to Ua m- 


4-dimensional covariant derivatives). The 3-D (spatial) covariant-derivative operators, on the 
other hand, will be D^ = (orthogonal to the Ua-field) and Da = (normal to Ua)- 

We should point out that the peculiar velocity field {va) lies in the 3-D space of the Ua-frame 
but not in that of tilted observers. In particular, one can easily show that UaV^' = ^ 0 and 

that hj’vh = Va + 'jv'^Ua 7 ^ Va- Of course, when <C 1 and 7 ci 1 , the tilt angle (/3) tends to 
zero and the aforementioned 3-D spaces essentially coincide (i.e. ea — Va - see Fig. [1]). 

3 Kinematics of peculiar motions 

Once the spacetime has been split into time and 3-D space, one can chose a family of observers 
and proceed to decompose every variable, every operator and every equation into their temporal 
and spatial parts relative to the chosen frame. For instance, the covariant derivative of the 
Ua-field splits into the irreducible components of its motion according to (e.g. see [lailH]) 

VbUa = ^Qhah +(Tah + VJah-AaUh- (3) 

Here, the (volume) scalar 0 = 'S/°'Ua describes the mean separation between the observers, with 
positive values of 0 indicating expansion and negative ones contraction. The symmetric and 
trace-free tensor aab = ^{h^a) represents the shear and the antisymmetric tensor ujab = Dj^Ua] 
the vorticity of the observers’ worldlines. Also, Aa = Ua = u^'^bUa is the 4-acceleration vector 
and vanishes when these worldlines are timelike geodesics. Similarly, the covariant derivative of 
the ha-field decomposes as 

Vfeha — ~ Qhab T dab “t“ iVab ^aVb j (4:) 
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where now 0 = V°‘Ua, ^ab = D(b'Wa)! ^ab = D[b'Wa] and Aa = u'a = u^^hUa- Starting from these 
definitions, one can derive the expressions relating the tilded variables to their non-tilded coun¬ 
terparts, in terms of the gradients of the peculiar velocity field. Treating the peculiar velocities 
as perturbations on an FRW background, we arrive at the linear transformation laws mm 

0 = 0 + = 0 + D“Ua (5) 

and 

Aa — Aa + Va T ~ QVa — Aa T Va T ^ 0^^a • ( 6 ) 

Note that, in deriving the above, one needs to remember that hab = hab + to first order 

in V. This guarantees that D“ua = D“ua at the same perturbative level. Also, using Eq. ([2]), one 
can easily verify that v'^ = Va to linear order in v. 

Relations ([5]) and ([ 6 ]) show that, even though the universe may be an exact FRW model 
relative to the tta-frame, it will not appear so to the tilted observers. The former, in particular, 
ensures that relatively moving observers measure different expansion rates. The latter, on the 
other hand, shows that we cannot simultaneously treat the worldlines of the relatively moving 
observers as timelike geodesics, even when the peculiar velocities are small. More specifically, 
in general we have Aa ^ 0 when Aa = 0 (or vice-versa), due to relative motion effects alone. 

4 Large-scale peculiar motions 

By construction, the expansion of the universe defines is a preferred reference frame, with respect 
to which one can define relative motions and measure peculiar velocities. This is the coordinate 
system of the smooth Hubble flow, namely the frame where the CMB dipole vanishes. In 
what follows we will always identify the Ua-field with the 4-velocity of an observer following the 
Hubble expansion. We will also assume that the aforementioned (fundamental) observer ‘sees’ 
the universe as an FRW cosmology, containing conventional pressure-free dust. This will be the 
‘background’ universe of the tilted observers, namely those with 4-velocity Ua, living in a typical 
galaxy (like our Milky Way) that moves relative to the CMB frame with peculiar velocity Va- 
The rest of this paper will look into the average kinematics of the aforementioned two families 
of observers. More specifically, we will focus on the rate of their mean (volume) expansion and 
on their acceleration/deceleration. Since we are dealing with non-relativistic peculiar motions, 
the volume scalars measured in the two frames are related by 018 ] 

© = 0 + , (7) 

with 'd = D“ua (see Eq. ([5]) above) H Note that 0 > 0 always due to the background expansion 
(i.e. that of the Ua-held). In addition, since we will be dealing with relatively slow peculiar flows, 
we will assume that I'd! <C 0. This makes Va and first-order perturbations. The above relation 
shows that 0 7 ^ 0 , that is the expansion rate of the tilted observer and that of the actual 
universe are generally different. In particular, when is positive, the tilted observers expand 

®By construction, the two volume scalars may be written as 0 = 3H and 0 = 3H, where H and H are the 
Hubble parameters measured in the Ua and the Ua frames respectively (e.g. see |12l I13 p 
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slightly faster than the actual universe, while in the opposite case (i.e. for < 0) their expansion 
is slightly slower. One therefore expects to see a difference in the acceleration/deceleration rates 
between the two frames. Indeed, taking the time derivative of the above, relative to the ha-frame, 
and keeping up to linear-order terms, we arrive at 


0' = 0 + h', 


( 8 ) 


where 0' = h“Va0, 0 = n“Va0 and h' = h“Vah [3 |8]. Accordingly, 0' ^ 0, with their 
actual difference depending on the sign and on the magnitude of hb Note that, although current 
observations strongly suggest that h <C 0 on scales greater than ~ 100 Mpc, we should not a 
priori extend this condition to their time derivatives. 

Different values for 0' and 0 imply that the deceleration parameters, as measured in the 
two frames, will differ as well. By dehnition, we have 


9 = - 



and 



(9) 


for the deceleration parameters in the tilted and the CMB frame respectively. Solving the above 
for 0' and 0 and substituting the resulting expressions into Eq. ([8]), gives 


1+q 


(1 + 9 ) 





( 10 ) 


This provides our hrst relation between the deceleration parameter measured in the rest frame 
of a typical galaxy (q) and that of the actual universe {q). Moreover, recalling that & = & + '& 
and performing a simple Taylor expansion reduces the above into 


1 + 9 


(1 + 9) 



02 ’ 
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after keeping up to n-order terms. Therefore, in the absence of peculiar motions (i.e. when 
= 0 = I?'), the two deceleration parameters coincide (as expected). Generally, however, q ^ q 
and their difference is mainly decided by i?', namely by the time evolution of 'd. This scalar, 
which has been dehned as i? = D“ua, describes the mean expansion/contraction of the peculiar 
flow0 

In principle, one could calculate the ratios 7?/0 and 'd'/Q"^ from the observations (recall that 
0 = 3H in the CMB frame) and then use expression (|lip to estimate the deceleration parameter 
in the rest-frame of the peculiar flow. In practice, however, the surveys only provide the mean 

^The irreducible kinematic variables of the bulk motion are obtained by decomposing the spatial gradient of 
the peculiar velocity field. In particular, relative to the tilted frame, we have 

DfjUa — 'dhab ^ab '^ab , (t2) 

with -d = D“wa, <fa6 = and S^ab ~ D[i,Wa]. Therefore, in analogy with 0 and 0, the scalar i? describes the 

average separation of the peculiar-flow lines, namely their expansion or contraction. Similarly, ^ab represents the 
‘pecnliar’ shear and ujab is the ‘peculiar’ vorticity relative to the tilted frame [13 US]. 
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velocity {{v)) of the bulk. One could then use the (approximate) relation |i?| ~ \d°‘Va\ — 3{v)/X, 
where A are the dimensions of the bulk motion, to obtain an estimate for the absolute value of 
the associated volume scalar. The problem is that, so far at least, the data cannot provide the 
sign of 7? on the scales of interest (i.e. for A > 100 Mpc). Even more difficult is to estimate the 
time derivative of (and therefore the ratio 'd 'directly from the observations. Therefore, 
in the next section, we will attempt to ‘ascertain’ i?' theoretically. 

5 The Raychaudhuri equation of the bulk flow 

The Raychaudhuri equation monitors changes in the mean separation between the worldlines of 
neighbouring observers, that is whether these expand/contract at an accelerated or decelerated 
pace. Although Raychaudhuri’s formula has been typically associated with timelike worldlines, 
it can be applied to curves of any type. For instance, changes in the average separation between 
the peculiar flow lines, which are spacelike curves, are monitored by what we may call the 
peculiar Raychaudhuri equation [16]. The latter can be obtained in the ‘traditional’ way, namely 
by applying the Ricci identities to the peculiar velocity field, or by taking the time derivative 
of 'd in the tilted frame. In either case, assuming non-relativistic peculiar flows in an FRW 
cosmology with dust, the calculation gives 

^' = -l0^ + bX, (13) 

to first order |16] . Note that in deriving the above we have used two sets of linear relations 
between the variables in two frames |T3] 

P = P, P = 0, qa = -pva and vfafe = 0. (14) 

Here, p is the density, p is the isotropic pressure, qa is the energy-flux and iTab is the anisotropic 
pressure of the matter, relative to the tilted observers, with p = 0 = qa = Tiab ia the CMB frame. 
Hence, the cosmic medium can be treated as pressure-free in both frames, to first approximation, 
but not as a perfect fluid. Indeed, in line with (114b ). the tilted observers ‘measure’ nonzero 
effective energy-flux due to their relative motion alone. Also, Eq. ([6]) gives 

^a = v'a + ^ QVa , (15) 

since Aa = 0 in the CMB frame. Therefore, Aa ^ 0 solely due to the observers peculiar flow. 

Essentially all the relative-motion effects discussed in the rest of this paper stem from the 
fact that Qa 7 ^ 0 and / 0 (see Eqs. (fTlh ) and ([Eb) respectively). The former shows that 
matter cannot be treated as a perfect fluid in the tilted frame and the latter implies that the 
worldlines of the tilted observers are no longer geodesics. Recall that = 0 = Aq in the 
CMB frame, relative to which the universe is an FRW cosmology filled with a single perfect 
fluid of zero pressure]! All these mean that linearising around a Friedmann background would 

®The fact that Aa = 0 in the CMB frame and that Aa 7^ 0 in the tilted frame, has immediate and important 
implications for the form of the associated Raychaudhuri equations. In particular, recalling that p = p and 
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generally produce different (linear) equations in the CMB frame and in the tilted frame, due to 
relative-motion effects alone (e.g. see footnote 5 and § [6] next). 

Looking at the first term on the right-hand side of Eq. m, we notice that the universal 
expansion always slows the bulk peculiar motion down. In particular, expanding local flows have 
> 0, in which case the sign of the aforementioned term is negative and therefore decelerates 
the peculiar expansion. Similarly, when dealing with contracting bulk motions (i.e. for -d < 0), 
the expansion term becomes positive and slows the peculiar contraction down. The role of the 
second term on the right-hand side of expression (|13p . however, is less straightforward to decode 
and in order to do so we will turn to cosmological perturbation theory. 


6 The role of the inhomogeneities 


Bulk peculiar flows are believed to be the result of the ongoing structure-formation process, 
which increases the local inhomogeneity and anisotropy of the universe. In fact, one can directly 
relate the peculiar velocities to inhomogeneities in the density of the matter. In the tilted 
frame, these are monitored by the dimensionless variable Aq = {a/p)DaP, which describes 
spatial variations in the density between neighbouring (tilted) observers (e.g. see |13]). Keeping 
up to u-order terms, gives Aq = {a/p)T)aP to linear order. Then, at the same approximation 
level, this variable evolves as m 

^'a = ~^a H-- + q -j (17) 

P \ ^ J P 


with ^ = aDa© describing inhomogeneities in the volume expansion with respect to the tilted 
frameO. Recalling that = —pva in our case (see (fUb l above), Eq. (flTl) transforms into 


a; = -Z, - a© 


+ 3®”“ 


+ aBa'd , 


(19) 


p = p = 0 to linear order (see Eqs. (I14h ') and (I14b i above), we have (in geometrised units) 

e = -ie'-ip and e' = -i02 _ ip + D“^, (le) 

relative to the Wa-frame and its tilded counterpart respectively. Given that that Aa = v^ + (0/3)va (see Eq. (IlSIl i. 
expression (llBb i shows that relative-motion effects, alone, can change the deceleration/acceleration rate of the 
tilted observers. This happens in the absence of pressure and constitutes a significant theoretical difference 
between our two-fluid study and the ‘standard’ single-fluid approaches. In fact, the aforementioned difference is 
the reason behind the relative-motion effects discussed in this paper (see §[7]and §l8] below). 

Note that an alternative way of obtaining expression m is to start from Eqs. (HSJ and then use the linear 
relations and m, while keeping in mind that <C 0. Finally, it is worth pointing out that combining 
Eqs. (11611 with definitions Q and then using ([7]), (1 1311 and (I15II leads to expression (Hill . 

®In the CMB frame, namely relative to the Ua-field, linear inhomogeneities in the density distribution of the 
(pressureless) matter evolve according to 

Aa = —Za , 

with Aa = oDaP and Za = aDa0 (e.g. see |13p. Note that, to first approximation, Aa yf Aa and Za yf Za 


( 18 ) 



since i? = T)°‘Va and p'va 
arrive at 


QpVa- Finally, taking the comoving divergence of the above, we 


A’ = -Z- 0^0 ( + - 0^ ) + , 


( 20 ) 


where A = aD“Aa and Z = aD“Za. Also, aD“A(j = A' to linear order. By construction, 
the scalar A monitors overdensities or underdensities in the matter distribution, relative to the 
tilted frame (e.g. see m)- Solving Eq. ([20]) for and substituting the resulting expression 
into the right-hand side of (fTHl) gives 






A' 

^ ©j ■ 


( 21 ) 


This result provides the linear evolution formula of r?, namely the linearised peculiar Raychaud- 
huri equation (see m for a more general expression), written in the frame of the tilted observers. 
Note the Laplacian term on the right-hand side of Eq. m, which introduces a scale-dependence 
that will prove crucial later. 

Our next step is to harmonically decompose the perturbations, by setting ^ = 'En'd{n)Q^^\ 
A = and Z = where Dai?(n) = 0 = DaA(„) = Da-^(n) and n represents 

the comoving eigenvalue of the harmonic mode. Also, are standard scalar harmonic func¬ 
tions, with = 0 and = —{n/a)‘^Q^^\ On using all of the above, expression ([2T]) 


takes the form 




/n\2 f ^ 

Va/ 0 a2 0 0 ^ ■ 


( 22 ) 


Finally, since An = ajn is the physical scale of the perturbation (i.e. the size of the bulk flow). 
An = 1/if = 3/0 is the Hubble radius and Ak = a/|if| (with A = ±1) is the curvature scale 
of the universe, we arrive at the expression 


d- 


(n) 


02 


2 

3 


1 





i {^ 

9 \Ak) 0 ^ 0 ^ ■ 


(23) 


which shows that the value of the dimensionless ratio is scale dependent. In fact, there 

are three scale-related thresholds in Eq. (1231) . corresponding to the size of the bulk flow (A^), to 
the Hubble radius {Ah) and to the curvature scale {Ak) of the universe. What is perhaps more 
important is that, due to the scale dependence of the above relation, small values of ??(„)/0 do 
not necessarily imply the same for the ratio ?9j^j/02, as seen from the tilted frame. When dealing 
with bulk flows much smaller than the horizon (i.e. those with Ah/A n 3> 1), in particular, the 
smallness of '(9(„)/0 does not guarantee the same for if(„)/0^- 


7 The scale-dependence of q 

By means of Eq. dm), the scale dependence of passes on to the deceleration parameter, 
as measured in the tilted frame. More specifically, substituting (I23h into the right-hand side of 
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([ni), gives 


1 + 9(n) 


(1 + g) 



+ 2 




+ i ^ ^ ^ 
0 3 \XkJ y e q 


(24) 


to first approximation. Before proceeding to extract some representative numerical estimates for 
q(n), let us first look at the last term on the right-hand side of the above. Today, observations 
strongly indicate that [Xh/XkY = |i- “ ^k\ ^ 10“^. In addition, <C 1 due to the 

anticipated large-scale homogeneity of the universe and the same could also be argued for the 
dimensionless ratio A|^^/0. Then, expression ([Mil reduces tclll 


1 + q(n) ^ (1 + ?) 1 - 2 - 




(n) 


0 


+ 2 


1 + 


KJ 




(n) 


0 


(25) 


0bservations of bulk peculiar motions typically cover scales of few hundred Mpc m 0], 
although there have been reports of peculiar flows on Gpc scales [6]. In all cases, the size of 
the drifting region lies inside the horizon (i.e. A„ < Xh) and the value of the peculiar velocity 
is a small fraction of the Hubble speed on the corresponding scale (i.e. 'd(„)/0 <C 1). On these 
wavelengths, expression (I25j) simplifies to 


1 

Qin) - ^ + 3 


Xh\^ "^{n) 
Xn) 0 


(26) 


providing a very simple relation between the deceleration parameter of the actual universe {q) 
and the one measured in the rest-frame of a typical galaxy (^(„)), which resides within a bulk 
flow of size Xn that moves relative to the smooth Hubble expansion. 

As expected, in the CMB frame where all peculiar velocities vanish identically, Eq. (j25p 
confirms that q = q on all scales, as expected. In the rest frame of the bulk motion, however, 
q ^ q and the difference also depends on the scale of the bulk flow. It comes to no surprise that 
the effect of the peculiar motion drops with increasing scale. In fact, on superhorizon lengths 
(with Xn/Xn 1) we find that q ^ q, which is something to be expected as well. The scale 
dependence seen in Eq. (1261) becomes the decisive factor as we move inside the Hubble radius and 
towards successively smaller lengths. In fact, on sufficiently small scales (i.e. for Xn/Xn ^ 1), 
even relatively slow bulk motions (i.e. those with /0 -C 1) can have a measurable effect on 
the deceleration parameter. Note, however, that one should not apply expression (l25]) to scales 
much smaller than 100 Mpc, since there the peculiar velocities are a considerable fraction of the 
Hubble speed (i.e. /0 ~ 1) and our approximation breaks down. 

The most intriguing feature of Eq. (j26p is that the relative motion might even change the 
sign of the deceleration parameter. This could happen when the peculiar flow is contracting 
(i.e. for -d^n) < 0) from now on, we will focus on contracting bulk flows only. Following our 

^Strictly speaking, the observations restrict the density contrast A(„) well bellow unity, but not its temporal 
derivative and therefore the dimensionless ratio A(^)/0. Nevertheless, given that (\h/\k)^ < 10“^, the product 
(A_ff/Ax)^(A(^)/0) should be negligible on all scales of interest. Put another way, (Air/AK)^(A(„)/0) 1, unless 

A(„)/0 ^ 10^ in regions of few hundred Mpc (or larger). The latter seems rather unlikely. 
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original assumption, that relative to the CMB frame the universe is an FRW model, we may set 
g = 1/2 in the right-hand side of Eq. (1261) . Also, since the background geometry is essentially 
Euclidean and the peculiar velocities are non-relativistic, we may write |7?| ~ d°'Va — 3{v)/X. 
Here, (v) is the mean magnitude of the peculiar velocity and A = is the size of the associated 
bulk motion (we have dropped the n-index for the economy of the presentation). Then, given 
that 0 = 3H, Eq. (1261) can be written as 


- 1 1 (v) 


(27) 


since we have assumed contracting peculiar flows (with 'i?/0 ~ —{v)/XH). Note also that 
q\ = in our new notation. Substituting into the right-hand side of the above the size and 
the mean velocity of a bulk flow, provides the value of the deceleration parameter measured by 
observers living within the flow. 


8 Estimating q 


In some cases the average peculiar velocity of the bulk flow may not be available. Also, sometimes 
the divergence of the observed mean bulk velocity may not accurately provide the volume scalar 
of the flow. It is conceivable, in particular, that |i?| < 3{v)/X. One could account for this 
by introducing a parameter (a), so that d = —3a{v)/X for contracting bulk flows, and then 
recasting expression (l27)l into 


1 1 /AhV a{v) 

2-n“J 


(28) 


with 0 < a < 1 in most cases. Next, we will employ the above and use peculiar-velocity 
measurements from recent surveys of large-scale bulk motions to obtain numerical estimates for 
the values of the deceleration parameter measured inside these flows. 

Eollowing (j28p . the value of qx depends stronger on the scale of the bulk motion than on its 
speed. In particular, the effects of relative motion on qx are more pronounced on scales well inside 
the Hubble radius. Such large-scale peculiar flows have been recently reported by several groups. 
Most of the surveys cover lengths close to lOO/Zi ~ 150 Mpc, though they have been claims of 
coherent flows with radii of Mpc-order (the so-called dark flows). The latter reports seem to be 
at odds with the data collected from the Planck mission (see however |10]1. but peculiar motions 
on scales smaller than 400 Mpc are not constrained by Planck. These surveys seem to agree on 
the direction of the bulk flow, but vary on the magnitude of the peculiar velocity. Some of the 
reported magnitudes are close to the general expectations of the ACDM scenario, though others 
are in excess of the predicted values. Here, we use representative peculiar-velocity magnitudes 
of large-scale bulk motions, to estimate the value of the deceleration parameter measured by 
observers living within the flow. In all cases we will assume that q = 1/2, 77 ~ 70 km/sec Mpc 
and Xh — 3 X 10^ Mpc. 

The most ‘conservative’ claims are probably those by Nusser and Davis, reporting mean 
peculiar velocities close to (v) ~ 250 km/sec on scales of roughly A ~ 150 Mpc [5]. This 
data, which are in agreement with the ACDM constraints, correspond to (Xh/X) ~ 20 and 
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Table 1: Representative values of the deceleration parameter (q), based on the peculiar velocities 
((u)) reported from several bulk-motion surveys, as measured by observers living inside these 
flows. Note that the ‘background’ universe is a decelerating Friedmann model (with q = 1/2) 
and in the observers’ rest-frame the value of the deceleration parameter changes entirely due to 
relative-motion effects. The numerical estimates follow from Eq. (j28p . for two different values 
of the a-parameter, assuming that Ff ~ 70 km/sec Mpc and Xh — 3 x 10^ Mpc. Note that the 
data from Nusser &: Davis (first line) are in agreement with the ACDM model. 


Survey 

A (Mpc) 

(v) (km/s) 

a = 1 

a = 1/2 

Nusser & Davis (2011) 

150 

250 

9x - -2-16 

q^ ~ —0.83 

Lavaux, et al (2010) 

150 

370 

q^ ~ -3.50 

q^ ~ -1.50 

Watkins, et al (2009) 

150 

400 

q^ ~ -4.83 

Qx - -2-16 

Colin, etal (2011) 

250 

260 

q^ ~ -0.22 

q^ ~ -1-0.14 

Abate (fe Feldman (2012) 

1000 

4000 

q^ ~ -hO.32 

q^ ~ -1-0.41 

Ade, et al (2014) 

>500 

<400 

Qx ^ +0.38 



{v)/XH ~ 0.02. Substituting these values into the right-hand side of Eq. (f28]l and assuming that 
a = 1, namely that |7?| = 3(u)/A, we find that q ~ —2.16 on the associated scale (see Table 1). 
When I'd! < 3{v)/X, the value of the local deceleration parameter increases. Setting a = 2/3 
and a = 1/2, for example, gives q = —1.27 and q = —0.83 respectively. On the same scales, 
Lavaux et al and Macaulay et al have reported larger peculiar velocities (around 370 km/sec) [3]. 
Using this survey, the deceleration parameter can go as low as g ~ —3.50. Much lower values 
for q are obtained from the surveys of Watkins et al [2], reporting peculiar velocities close to 
400 km/sec on scales of approximately 150 Mpc (see Table 1). 

We should emphasise again that none of the aforementioned values of the deceleration pa¬ 
rameter are ‘real’, but they are simply the result of the observers’ relative motion with respect to 
the smooth Hubble expansion. In the CMB frame, the peculiar velocities vanish identically and 
q = q = 1/2, as expected in an FRW universe. Intuitively speaking, observers residing inside a 
contracting bulk flow could be misled to believe that the background expansion accelerates, in 
the same way passengers in a vehicle that slows down might think that the overtaking cars are 
speeding up. 

Following (j28h , the impact of the relative motion on q weakens significantly as we move on to 
progressively larger wavelengths, where the ratio {Xh/X)‘^ drops sharply. Colin et al, for example 
have reported peculiar velocities around 260 km/sec on scales close to 250 Mpc [3j. Substituting 
the findings of this survey into the right-hand side of Eq. (j28p . we find that q can take only 
marginally negative values. Also, on Gpc-order scales, Kashlinsky et al and Abate et al have 
reported bulk motions as fast as 1000 km/sec and 4000 km/sec respectively [6]. Nevertheless, 
despite the large amplitude of the peculiar velocities, their effect on the value of the deceleration 
parameter, as measured within these flows, is relatively weak. Using the velocities reported by 
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Abate et al, in particular, expression (l28t) gives q ~ +0.32 on the corresponding lengths (for 
a = 1 - see Table 1). 

Finally, the analysis of the Planck data seems to set an upper limit of {v) ~ 400 km/sec 
on the mean velocity of bulk motions larger than 350//i ~ 500 Mpc [9]. Using this constraint, 
Eq. (1281) puts a lower limit of g ~ +0.38 on the corresponding scale (see last line in Table 1). 

9 Discussion 

Observers moving relative to each other could have a different view of ‘reality’, and they might 
also interpret their observations differently. In this study we have assumed two families of 
observers with a non-relativistic peculiar velocity between them. The first family was identified 
with the smooth expansion of a dust-dominated Friedmann universe and the second (the so- 
called ‘tilted’ observers) with a typical galaxy moving relative to the Hubble flow. Our aim 
was to investigate whether relative motion modifies the kinematics of the drifting observers. 
More specifically, we wanted to look for differences in the value of the deceleration parameter, 
as measured in the aforementioned two frames, that were solely due to relative-motion effects. 
We were primarily motivated by resent surveys, reporting coherent bulk motions with sizes of 
few hundred Mpc and peculiar velocities in excess of those anticipated by the ACDM paradigm. 
In agreement with mi, the deceleration parameter in the rest-frame of large-scale bulk flows 
was generally found to differ from that of the background universe. Extending that work, we 
realised that the aforementioned difference also depends on the scale of the bulk motion and 
not only on the magnitude of the peculiar velocity. In particular, the smaller the size of the 
drifting domain (relative to the Hubble radius), the stronger the relative-motion effects on the 
local deceleration parameter. Moreover, observers inside slightly contracting bulk flows could 
measure a negative deceleration parameter in their own rest-frame, while the universe they live 
in is actually decelerating. 

None of the aforementioned effects is real of course, given that they all vanish in the CMB 
frame of the smooth Hubble expansion. Just like the CMB dipole, all of the above are nothing 
else but mere artifacts of the observers’ relative motion. Nevertheless, when not properly taken 
into account, such relative-motion effects could potentially lead to the misinterpretation of the 
observations. To avoid this, one needs to abandon the single-fluid models and allow for a second 
set of observers moving relative to the Hubble flow. It is then straightforward to show that to 
these observers the cosmic medium no longer looks like a perfect fluid, even when their peculiar 
velocity is small. As a result, the Raychaudhuri equation in the rest-frame of a large-scale 
bulk motion has an additional term, which carries the relative-motion effects (see Eq. (|16b l 
in footnote 4). The latter are immediately transferred to the local deceleration parameter (see 
Eqs. ()24p - (l28p l. making it different from that of the actual universe (sometimes even in the sign). 

Following Table 1, the deceleration parameter measured inside slightly contracting bulk 
motions of few hundred Mpc is particularly voulnerable to relative motion effects. Adopting the 
peculiar velocities claimed by Watkins et al, which are in tension with the ACDM paradigm, 
we obtained values as low as g ~ —4.83 on scales around 150 Mpc. On the same wavelengths, 
the more conservative bulk velocities reported by Nusser et al, which are within the ACDM 
limits, gave a lowest value of g ~ —2.16. Analogous numerical estimates for the deceleration 
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parameter were claimed in m- There, the authors used a phenomenological two-parameter 
function for their deceleration parameter to reconstruct its evolution over the last couple of 
billion years, assuming a dark-energy dominated (exact) FRW universe,. Here, we have arrived 
at a qualitatively similar result, working in a tilted (perturbed) Friedmannian cosmology with 
pressure-free matter. Finally, according to (|28p . the effects of relative motion on the deceleration 
parameter inside the bulk flow, drop sharply with increasing scale. This means that q quickly 
approaches its ‘background’ value (i.e. q ^ q = 1/2), as we move on to progressively longer 
wavelengths. To the unsuspecting observer living inside the bulk flow, however, this may give 
the false impression that the universal acceleration is a relatively recent event. 

Is it then possible that the recent accelerated expansion of the universe could be a mere 
illusion and a relative-motion artifact? Although this may not be inconceivable, in principle, it 
is also rather premature to make such a claim. What we could say at this point is that there are 
many attractive aspects in this alternative, assuming that it works. There is no need for dark 
energy, it is not necessary to modify gravity and there is no reason to abandon the Friedmann 
models. Everything happens within standard physics and conventional cosmology. In addition, 
since on the scales of interest the observed peculiar velocities are much smaller than the Hubble 
speed, we have not yet included any nonlinear effects. There seems to be a key requirement 
however. In order to create the (false) impression of accelerated expansion, the bulk flow we live 
in must be slightly contracting (on a background that is expanding). In a sense, observers living 
inside a contracting bulk flow may think that the background expansion is accelerating, just 
like passengers in a car that slows down are sometimes mislead to believe that the overtaking 
vehicles are speeding up. At the moment, we cannot measure volume changes (i.e. the value and 
the sign of d) for bulk motions of 100 Mpc-order or greater. Nevertheless, provided there is no 
bias in favour of contraction or expansion on these wavelengths, the chances of living within a 
slightly contracting large-scale peculiar flow should be around 50%. 

There are caveats as well of course. For example, we have assumed that the volume scalar 
(i?) of the bulk flow can be obtained from the its observed mean velocity ((u)), by means of 
the approximate formula 'll ~ 3a{v)/X, with 0 < a < 1. It is conceivable, however, that the 
a-parameter might take very small values. For instance, if the bulk flow is incompressible, then 
a, 1 ? —0 and the relative-motion effects on q become negligible for all practical purposes. A 
more subtle caveat lies perhaps in our assumption that the ratios A'^^^/0 and Z(„)/0 inside the 
last term of Eq. <\2^ are much smaller than unity and therefore their effect on q is irrelevant. It 
is not a priori inconceivable, however, that in the rest-frame of the bulk motion these ratios may 
not be small and that their effect on q may not be negligible. In particular, the density contrast 
A could vary rapidly with time in the tilted frame, in which case the ratio A'/0 may take 
relatively large values. Having said that, the coefficient {Xh/Xk)"^ of the aforementioned last 
term is very small. More specifically, {Xh/Xk)‘^ ^ 10“^, according to the current observational 
constrains. Therefore, unless A'^^^/0 -|- > 10^ in the tilted frame, the last term on the 

right-hand side of Eq. (j24|] will be irrelevant. It should also be noted that, if the inhomogeneous 
(i.e. the last) term dominates the right-hand side of (|24p . the value of the deceleration parameter 
in the tilted frame will still differ from that of the actual universe (i.e. q ^ q)- In that case, 
however, the conditions leading to negative q may not be the same with the ones discussed in 
§ [3 and § [8] earlier. 0verall, the only way of having g ~ g, is when the last two terms on the 
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right-hand side of Eq. (1241) essentially cancel each other out. It is conceivable that some of the 
these issues may be resolved when nonlinear effects are included into the study. 

We would like to close our discussion by pointing out that we have assumed observers with 
individual peculiar velocities equal to the average velocity of the bulk flow. In general, however, 
a typical observer inside the bulk will have peculiar velocity close (in terms of both magnitude 
and direction) to the mean, but not equal to it. Put another way, typical observers will have their 
own (small) peculiar velocities relative to the average motion of the bulk flow. To these observers, 
the distribution of q will not be entirely isotropic, but it should demonstrate a weak dipole-like 
anisotropy analogous to that seen in the CMB. Recall that the latter is caused by the motion of 
the bulk with respect to the smooth universal expansion. In a similar way, the aforementioned 
g-anisotropy will be due to the observer’s peculiar flow relative to the mean motion of the bulk. 
Assuming that the difference between the two velocities is small, the anisotropy in the observed 
distribution of q will be weak and the associated dipole axis should lie close to that of the CMB 
(see [8] for further discussion). Interestingly, there have been surveys arguing that such a weak 
dipolar anisotropy might actually exist in the supernovae data. In other words, there have been 
claims suggesting that our universe is accelerating slightly faster in one direction and equally 
slower along the opposite (e.g. see [H]). Moreover, the dipole axis seems to lie fairly close to 
that of the CMB. If these reports are confirmed, one might want to take a closer look at the 
scenario outlined here as an alternative explanation of the recent universal acceleration. So far, 
however, the presence of a weak dipolar axis in the supernovae data remains ambiguous. 
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